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SECTION  I 
INTRODUCTION 


BACKGROUND  AND  SUMMARY 


Theoretically  derived  shapes  of  collision-broadened  H2O  lines  have 
been  well  verified  in  regions  within  a  few  cin"l  of  the  line  centers. 
However,  it  is  well  knownl”^  that  the  absorption  by  the  extreme  wings 
of  the  lines  cannot  be  accurately  calculated  solely  on  the  basis  of  the 
known  line  parameters:  position,  intensity  and  width.  The  atmospheric 
windows  centered  near  1000  cm"l  and  2600  cm~l  are  two  examples  of 
spectral  regions  where  HpO  absorption  is  quite  different  from  that 
calculated  from  theoretical  shapes.  A  major  reason  for  the  disagreements 
with  the  absorption  theory  is  the  lack  of  understanding  of  the  shapes  of 
the  extreme  wings  of  the  absorption  lines.  Measurements  of  the  weak 
absorption  in  the  major  windows  show  that  the  extreme  wings  of  the  lines 
deviate  significantly  from  the  theoretical  shapes  beyond  approximately  50 
cm"^  from  the  line  centers. 


In  a  previous  report^,  we  have  shown  evidence  that  the  actual 
shapes  of  the  H2O  lines  also  deviate  from  theoretical  at  closer 
distances  from  the  line  centers,  probably  as  close  as  5  or  10  cm~l. 
Most  workers  in  this  field  had  previously  assumed  that  the  theoretical 
shapes  were  still  correct  at  greater  distances  than  this  from  the 
centers.  The  previous  report^  that  indicates  deviation  from 
theoretical  shapes  in  the  intermediate  wings  deals  with  absorption  in 
several  narrow  windows  between  the  strong  lines  centered  from 
approximately  1300  cm'^  to  2000  crn'l.  The  present  report  summarizes 
a  similar  study  of  approximately  50  narrow  windows  between  strong  lines 
from  3000  cm'l  to  4200  cm”l. 


The  general  results  of  the  present  work  are  consistent  with  those  of 
the  similar  study  in  a  different  vibration-rotation  band.  In  all  of  the 
narrow  windows  within  the  main  parts  of  the  bands,  the  observed 
absorption  is  greater  than  that  predicted  from  the  Lorentz  line  shape  and 
the  AFGL^  parameters.  The  relative  difference  between  experimental  and 
theoretical  results  is  greater  for  self  broadening  than  for  N2 
broadening.  This  relative  difference  is  also  greater  in  regions  where 
most  of  the  absorption  is  due  to  lines  centered  more  than  10  cm"^  away 
than  it  is  in  regions  very  near  the  strong  lines. 


1.  D.  E.  Burch,  AFGL-TR-81-0300,  ADA112264,  Final  Report,  AFGL  Contract 
No.  FI 9628- 7 9- 0041  (1982). 

2.  D.  E.  Burch  and  D.  A.  Gryvnak,  Continuum  Absorption  by  H2O  Vapor  in 
Atmospheric  Water  Vapor  (A.  Deepak,  T.  D.  Wilkerson,  and  L.  H. 
Ruhnke,  eds).  Academic  Press,  New  York  (1980). 

3.  D.  E.  Burch,  SPIE  Proceedings  277,  28  (1981). 

4.  R.  A,  McClatchey.  W,  S.  Benedict,  S.  A.  Clough,  D.  E.  Burch,  R.  F. 
Calfee,  K.  Fox,  L.  S.  Rothman,  and  J.  S.  Garing,  AFCRL  Atmospheric 
Absorption  Line  Parameters  Compilation,  AFCRL-TR-73-0096,  U.  S.  Air 
Force  (1973).  (Avalable  from  NTIS),  AD762904. 
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We  attribute  the  "excess"  absorption  at  the  wavenumbers  of  minimum 
absorption  to  a  super-Lorentzian  line  shape  beyond  about  5  or  10  cm"l 
from  the  line  centers.  It  cannot  be  due  entirely  to  errors  in  the 
assumed  intensities  or  widths  of  the  lines.  Errors  of  this  type  would 
cause  deviations  between  the  theoretical  and  experimental  results  at 
positions  close  to  the  strong  lines.  Such  errors  were  not  observed. 

The  excess  absorption,  which  we  call  empirical  continuum,  does  not 
exhibit  rapid  variations  with  changing  wavenumber  throughout  the  band, 
but  it  is  greatest  in  regions  that  contain  the  strongest  lines.  It  is 
too  small  to  be  checked  very  near  the  strong  lines,  where  the  absorption 
is  dominated  by  the  lines  themselves.  Thus,  ignoring  the  empirical 
continuum  in  a  line-by-line  calculation  leads  to  significant  errors  only 
in  the  windows  a  few  cm"^  from  the  lines  producing  most  of  the 
absorption.  The  average  transmittance  over  wide  spectral  intervals  is 
not  strongly  dependent  upon  the  empirical  continuum  when  the  average 
transmittance  is  high.  But  the  errors  in  calculated  average  transmittance 
may  be  large  for  long  paths  over  which  the  only  significant  transmittance 
is  in  the  windows  where  the  empirical  continuum  plays  a  major  role. 

DEFINITIONS,  SYMBOLS  AND  PARAMETERS 

At  pressures  corresponding  to  the  lower  atmosphere,  the 
vibration-rotation  lines  of  infrared  H2O  vapor  are  collision 
broadened.  It  has  been  well  established  that  the  shape  of  one  of  these 
lines  within  a  few  cm"l  of  its  center  is  given  by  the  simple  Lorentz 
equation. 

k  =  S  OL  ,  (simple  Lorentz)  (1) 

where  k  is  the  absorption  coefficient  at  wavenumber  v,  and  is  the  line 
center.  It  has  also  been  well  established  that  the  line  intensity  S  is 
essentially  independent  of  the  half  width  a  ,  which  is  proportional  to 
pressure.  For  a  line  shape  to  be  exact,  the  line  intensity  S  must  be 
related  to  k  by 

S  =  j'J  d  V  .  (2) 

The  simple  Lorentz  shape  given  by  Equation  (1)  is  known  to  be  invalid 
in  the  far  infrared  and  millimeter  regions,  or  in  the  extreme  wings  of 
any  lines  where  |v-Vo|  is  no  longer  much  less  than  a  ,  References  5  and 
6  discuss  more  appropriate  theoretical  shapes  for  all  spectral  regions 
and  their  deviations  from  experimental  results.  No  line  shape  derived 
strictly  on  theoretical  considerations  agrees  with  experimental  results 
in  all  of  the  spectral  regions  from  the  visible  to  the  microwaves. 


5.  S.  A.  Clough,  F.  X.  Kneizys,  R.  Davies,  R.  Gamache,  and  R.  Tipping, 
Theoretical  Line  Shape  for  H2O  Vapor:  Application  to  Continuum,  in 
Atmospheric  Water  Vapor  (A.  Deepak,  T.  D.  Wilkerson,  and  L.  H. 
Ruhnke,  eds.),  Academic  Press,  New  York  (1980). 

6.  S.  A.  Clough,  F.  X.  Kneizys,  L.  S.  Rothman,  and  W.  0.  Gallery,  SPIE 
Proceedings  272  (1981). 
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This  report  presents  new  data  on  H2O  vapor  absorption  between 
approximately  3000  cm"l  and  4200  cm“^,  and  compares  the  results  with 

a  previous  similar  study  of  the  1250-2200  cm"i  region.  Essentially 

none  of  the  absorption  in  either  of  these  intervals  is  due  to  lines 

centered  outside  of  the  interval.  In  calculating  the  absorption  within 
this  region,  the  results  obtained  by  summing  the  contributions  by  lines 
with  the  Lorentz  shape  are  very  nearly  the  same  as  those  obtained  by 
using  the  more  complex  form  developed  in  References  5  and  6.  Therefore, 
for  simplicity,  we  have  chosen  to  use  the  more  common,  simple  Lorentz 
form  for  calculations  to  compare  with  experimental  results. 

The  emphasis  in  this  report  is  on  the  absorption  in  approximately  50 
very  narrow  windows  between  the  lines.  At  the  wavenumbers  selected  for 
study,  very  little  of  the  absorption  is  due  to  an  absorption  line 
centered  closer  than  1  cm“l.  In  every  case,  the  measured  absorption  is 
greater  than  that  predicted  by  summing  the  calculated  contributions  by 
the  lines.  The  excess  absorption  is  the  main  subject  of  this  study. 

The  absorption  coefficient  for  a  single  line  is  related  to  the  true 
transmittance  T'  of  the  line  that  would  be  observed  with  infinite 
resolving  power  (zero  slit  width)  by: 


k  =  (-1/u)  In  T' 

The  absorber  thickness  u  is  given  by: 
u(molecules  cm“2)  *  p(atm)  7.34  x  10^^  L(cm)/9 


(3) 


(4) 


where  0  is  the  temperature  in  Kelvin.  L  is  the  geometrical  path  length 
of  the  radiation  through  the  sample  consisting  of  H^O  \apor  at  partial 
pressure  p.  It  follows  that  the  units  for  S,  defined  by  Equation  (2), 
are  molecules"^  cm^  cm’l. 

The  absorption  coefficient  given  by  any  theoretical  shape  for 
collision-broadened  lines  is  proportional  to  a  in  the  extreme  wings 
where  v-vq  >  >  a  .  This  proportionality  relationship  appears  experi¬ 
mentally  to  be  valid,  even  when  the  coefficient  is  greatly  different  from 
that  given  by  any  of  the  theoretical  expressions.  The  halfwidth  of  a 
collision-broadened  line  for  a  mixture  of  H^O  in  at  a  fixed 
temperature  is  given  by: 


0  _ 
a  S  P 


■a  N  PN 


(5) 


The  partial 
represented  by  p 


in 


are 


atmospheres  of  the  H^O  and  N2 
p|.,  respectively.  N2  is  used  to  simulate  dry 
^  are  assumed  to  be  proportional  to  as 


pressures 
and 

air.  Both  05  and  q  ^ 
recommended  by  Reference  4  for  H2O  lines.  The  equivalent  pressure  P 
is  proportional  to  regardless  of  the  relative  values  of  p  and  pfj. 


Pe  = 


Bp  +  PN  =  (B-1)  p  +  P 


(6) 


•  •  -A 


.“Ai' 


V*. 


where  P  is  the  total  pressure,  and 
B  =  a 
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(7) 


B  for  H2O  has  been  determined  experimentally  to  be  approximately  5; 
this  value  is  used  in  all  of  our  calculations  reported  here. 

The  total  absorption  coefficient  K  at  a  given  wavenumber  due  to  many 
lines  is 

K  =  5  ki  =  (-1/u)  In  T',  (8) 

where  each  represents  the  coefficient  due  to  an  individual  line,  u 
is  given  by  Equation  (4),  and  T'  is  the  true  transmittance  that  would  be 
observed  with  infinite  resolving  power  (zero  spectral  slit  width).  The 
present  work  deals  with  collision-broadened  lines,  each  of  v/hich  can 

conveniently  be  treated  as  having  two  regions,  one  near  the  center  and 
one  in  the  wings.  Thus,  we  write  K,  the  total  absorption  coefficient,  as 

K  =  -d/u)  In  T'  =  K(local)  +  K(wing)  (9) 

At  any  wavenumber,  K(local)  represents  the  contribution  due  to  nearby 
lines  for  which  [v-Vg  |  is  not  large  compared  to  a.  K(wing)  is  then 
defined  as  the  contribution  due  to  the  wings  of  lines  for  which  |v-vo|>>cf 
The  present  effort  deals  with  selected  wavenumbers  at  which  K( local)  is 
much  less  than  K(wing).  This  makes  it  possible  to  make  quantitative 

comparisons  between  the  experimental  and  calculated  values  of  K(wing). 
As  discussed  above,  the  main  purpose  of  this  study  is  to  compare  the 

experimentally  observed  wing  absorption  with  that  calculated  by  using  the 
simple  Lorentz  line  shape  given  by  Equation  (1).  The  K(local)  portion  of 
the  absorption  is  small,  but  significant,  and  is  assumed  to  be  given 
accurately  by  the  simple  Lorentz  line  shape.  The  method  of  accounting 
for  the  K(local)  contribution  is  explained  in  Section  3. 

All  of  the  experimental  values  of  window  absorption  measured  during 
the  present  study  are  greater  than  the  corresponding  calculated  values. 
We  define  the  difference,  or  "extra  absorption",  as  an  empirical 
continuum: 

ec  =  K(exp)  -  K(calc).  (10) 

Because  of  the  assumption  that  the  K(local)  contributions  to  the 
absorption  are  accounted  for,  the  empirical  continuum  is  equivalent  to 

the  differences  in  wing  absorption:  K(wing,  exp)  -  K(wing,  calc).  The 

results  indicate  that  the  extra  absorption  as  defined  by  the  above 
equations  does  not  vary  rapidly  over  spectral  intervals  a  few  cm‘1 
wide;  therefore,  the  word  continuum  is  appropriate. 

In  a  sample  of  pure  H2O  vapor  with  no  other  gas  present,  all  of  the 
line  broadening  is  due  to  collisions  of  the  absorbing  H2O  molecules 
with  other  H2O  molecules.  Such  H2O-H2O  collisions  lead  to  self 
broadening.  In  samples  of  H2O  +  N2»  both  self  broadening  and  N2 
broadening  of  the  H2O  lines  take  place,  and  each  contributes  to  the  K 
(wing)  coefficient. 

From  Equation  (1)  we  see  that  k,  the  coefficient  due  to  a  single 
line,  is  proportional  to  a  in  the  wings  where  |v-Voj>>a.  Thus,  K(wing), 
which  is  due  to  the  contributions  by  the  wings  of  all  the  lines,  is  made 
up  of  two  components,  one  proportional  to  p  and  one  proportional  to 
PN.  We  define  two  normalized  empirical  continuum  coefficients  eC9 


and  ®Cjj  that  represent  self  broadening  and  N2  broadening, 

respectively.  The  superscript  0  represents  a  normalization  to  1  atm 
pressure  of  H2O  vapor  or  N2.  Subscripts  s  and  N  denote  self 
broadening  and  N2  broadening,  respectively. 


For  self  broadening  only. 


p(atm) 


Ks(exp)  -  Ks(calc) 


(11) 


K5(calc)  is  determined  by  summing  the  contributions  of  all  of  the  lines 
with  the  following  assumptions:  that  the  line  shapes  are  given  by 


Equation  (1);  S,  Vq  ,  and  ot^ 
parameters  compilation(^) ;  and 
Equations  5  and  7. 


are  the  values  given  by  the  AF6L  line 
=  Son  p(atm),  in  accordance  with 


a: 


o  - 


For  a  mixture  of  H2O  +  N2. 
ec  =  eCs+eCN  =  P  ®C°3+ 


(12) 


where  =  Kf((exp)  -  Kfj(calc).  (13) 

KN(calc)  is  determined  by  summing  the  contributions  of  all  of  the  lines 
by  the  same  method  used  to  calculate  Ks(calc),  except  that  the  value 
used  for  half  width  is  aN. 

The  above  definitions  of  transmittance  T',  k,  K,  and  ®C  correspond 
to  a  single  wavenumber;  i.e.,  to  measurements  that  would  be  made  with 
infinite  resolving  power.  The  spectrometer  used  in  the  measurements 
passes  a  finite  spectral  interval,  typically  between  0.5  and  1  cm"l; 
thus  the  observed  transmittance  T(exp)  differs  from  the  true 
transmittance  T'.  The  method  described  in  Section  3  to  account  for  the 
finite  spectral  slitwidth  involves  calculations  of  T'(calc)  and  T(calc). 
At  a  fixed  wavenumber,  T'(calc)  for  a  given  value  of  u  is  calculated  by 
use  of  Equation  (8).  Values  of  T'(calc)  are  calculated  at  many 

wavenumbers  with  the  intervals  between  points  of  calculation  small  enough 
that  the  spectral  structure  is  retained.  These  values  of  T'(calc)  are 
then  convolved  with  a  slit  function  similar  to  the  known  experimental 
slit  function  to  determine  values  of  T(calc).  These  values  of  T(calc) 
are  compared  to  values  of  T(exp)  to  determine  the  deviation  of  the 

theoretical  results  from  the  experimental  results.  This  procedure  is 
illustrated  with  typical  spectra  in  the  discussion  of  data  reduction  in 
Section  3. 


After  determining  the  continuum  coefficient  at  each  wavenumber  studied, 
we  used  its  value  to  find  the  corrected  value  for  the  total  absorption  co¬ 
efficient.  We  then  compared  these  corrected  values  to  corresponding  values 
calculated  by  Air  Force  Geophysics  Laboratory  scientists  using  the  well 
known  FASCOOE®  computer  program. 
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SECTION  2 


EXPERIMENTAL 

The  experimental  methods  are  similar  to  those  we  have  used  for  several 
investigations  of  absorption  by  H2O  and  other  atmospheric  gases.  Gas 
samples  are  contained  in  either  of  two  multiple-pass  absorption  cells  that 
provide  optical  paths  from  approximately  2  meters  to  1000  meters. 
Spectral  transmittances  are  measured  with  a  grating  spectrometer  mounted 
in  a  tank  that  can  be  evacuated  to  avoid  absorption  by  atmospheric  gases 
in  the  optical  path  outside  of  the  sample  cell. 

SAMPLING 

All  of  the  new  data  reported  here  have  been  obtained  with  samples  of 
either  pure  H2O  vapor  or  of  H2O  +  N2  at  room  temperature,  296K.  The 
typical  sampling  procedures  are  as  follows; 

The  cell  is  evacuated  to  less  than  0.1  torr  immediately  before  H2O 
vapor  is  added.  Distilled  liquid  water  is  contained  in  a  small  "boiler" 
attached  to  the  sample  cell  by  a  short  piece  of  tubing  with  a  valve.  The 
tubing  is  connected  to  the  liquid  container  above  the  liquid  level  so  that 
only  the  vapor  can  pass  through  the  tubing  to  the  cell.  After  the  valve 
is  opened  to  allow  H2O  vapor  to  enter  the  cell,  the  liquid  is  heated  to 
speed  up  the  evaporation  and  the  transfer  of  vapor  from  the  container  to 
the  sample  cell.  Approximately  10-15  minutes  is  normally  required  to  add 
H2O  to  a  pressure  of  15  torr  (approximately  0.02  atm).  The  pressure  is 
allowed  to  stabilize  before  spectral  measurements  are  made. 

Another  tube  with  a  valve  connects  the  upper  portion  of  the  liquid 
container  to  a  vacuum  pump  that  removes  any  air  from  the  container  before 
the  valve  is  opened  in  the  tubing  to  the  sample  cell.  The  tubing  to  the 
pump  and  sample  cell  are  electrically  heated  to  avoid  condensation  in 
them. 

Sample  pressures  below  approximately  0.1  atm  are  measured  with  an  oil 
manometer;  pressures  from  0.1  atm  to  2  atm  are  measured  with  a  mercury 
manometer.  A  Dubrovin  gauge  serves  for  pressures  between  2  atm  and  15 
atm.  The  manometers  and  gauge  are  at  room  temperature,  making  it 
necessary  to  use  extra  care  to  avoid  condensation  in  them  while  studying 
samples  at  pressures  above  the  H2O  vapor  pressure  at  room  temperature. 
Normally,  the  manometer  or  gauge  and  the  tubing  connecting  them  to  the 
sample  cell  are  filled  with  dry  N2  to  a  pressure  slightly  above  the 
sample  pressure;  the  valve  is  then  opened,  allowing  a  small  amount  of  dry 
N2  to  flow  into  the  cell  until  the  pressure  comes  to  equilibrium. 
Because  of  the  very  large  volume  of  the  sample  cell,  the  small  amount  of 
dry  N2  entering  from  the  manometer  line  makes  a  trivial  change  in  the 
sample  pressure. 

Samples  of  H2O  +  N2  are  made  by  first  adding  the  H2O  to  the 
evacuated  cell  and  allowing  the  pressure  to  come  to  equilibrium.  If  the 
H2O  is  added  and  allowed  to  remain,  the  pressure  decreases  by  20-30% 
during  the  first  24  hours  and  may  continue  to  decrease  slowly  for  another 
I  or  2  days  because  of  the  slow  absorption  of  H2O  on  the  walls  of  the 
cell.  The  H2O  pressure  is  normally  brought  to  near  equilibrium  much 


10 


more  quickly  by  first  adding  H2O  to  a  pressure  10-20%  higher  than  the 
desired  equilibrium  pressure.  The  gas  is  allowed  to  remain  for  several 
minutes,  then  part  of  it  is  pumped  out  until  the  pressure  is  approximately 
equal  to  the  desired  equilibrium  pressure.  If  this  sample  is  allowed  to 
remain,  the  pressure  does  not  normally  change  more  than  a  few  percent 
during  the  next  several  hours.  The  pressure  of  a  pure  H2O  sample  is 
measured  again  at  the  time  of  a  transmittance  measurement. 

Dry  N2  from  commercial  cylinders  is  slowly  added  to  the  H2O  in  the 
sample  cell  through  4  different  ports  distributed  evenly  along  the  length 
of  the  cell.  While  the  N2  is  entering  the  cell,  mixing  fans  inside  the 
cell  are  turned  on  periodically  to  ensure  good  mixing.  Immediately  before 
the  spectral  data  for  the  H2O  +  N2  sample  are  obtained,  the  H2O 
partial  pressure  is  measured  with  a  dewpoint  meter  inside  the  sample 
cell.  The  results  are  compared  with  the  H2O  pressure  measurements  made 
before  adding  the  N2:  one  with  the  dewpoint  meter  and  one  with  the 
manometer.  The  dewpoint  is  measured  periodically  while  a  sample  of  H2O 
+  N2  ''S  ^0  monitor  the  H2O  partial  pressure.  Any  small 
changes  that  occur  are  accounted  for  during  the  data  reduction. 

Possible  errors  in  transmission  measurements  due  to  adsorption  of 
water  on  the  optical  surfaces  in  the  longer  cell  are  greatly  reduced  by 
heating  the  windows  and  the  internal  mirrors  to  approximately  3  K  above 
the  gas  temperature.  Transmission  measurements  made  with  pure 
non-absorbing  N2  indicate  that  there  is  negligible  effect  due  to  any 
turbulence  or  temperature  non-uniformities  caused  by  the  heated  optical 
components. 

TRANSMISSION  MEASUREMENTS 

Infrared  radiation  from  a  Nernst  glower  is  chopped  and  directed 
through  a  window  into  the  multiple-pass  absorption  cell  that  contains  the 
sample  to  be  studied.  After  the  beam  of  energy  has  passed  back  and  forth 
the  desired  number  of  times  in  the  cell,  it  is  directed  out  through  an 
exit  window  to  a  grating  spectrometer. 

The  spectrometer  is  used  either  to  scan  a  portion  of  the  spectrum  of 
interest  or  to  remain  set  so  it  will  pass  only  a  narrow  spectral  interval 
at  one  of  the  wavenumbers  where  the  continuum  absorption  is  being 
measured.  These  narrow  intervals  are  chosen  between  H2O  absorption 
lines  where  there  is  little  interference  by  nearby  lines  in  the  continuum 
measurement. 

The  transmittance  of  a  sample  is  determined  by  comparing  the  level  of 
the  infrared  signal  measured  with  the  sample  in  the  cell  to  that  measured 
when  the  cell  is  evacuated.  The  background  spectrum,  or  background  signal 
leve’  at  a  few  specific  wavenumbers,  is  typically  obtained  before  and 
after  the  measurement  on  the  gas  sample.  The  number  of  passes  of  the 
sample  cell  can  be  changed  without  disturbing  the  sample,  making  it 
possible  to  investigate  samples  of  different  absorber  uhickness  from  a 
single  batch  of  gas.  When  this  procedure  is  followed,  background  data  are 
obtained  with  the  evacuated  cell  at  all  of  the  same  path  lengths  used  for 
the  samples. 


SECTION  3 


RESULTS  AND  DISCUSSION 

Experimental  data  presented  in  this  section  for  the  3000-4200  cm“l 
region  are  from  two  sources.  The  data  on  self  broadening  have  been 
obtained  recently  under  the  present  contract,  whereas  those  data  on  N2 
broadening  are  from  one  of  our  previous  reports^.  The  analysis  of  the 
data  is  similar  to  that  used  previously  by  us^'^  for  data  obtained  in 
narrow  windows  of  the  strong  H2O  band  between  1250-2200  cm"^.  Some  of 
the  results  of  the  previous  study  are  repeated  in  this  section  for 
comparison  with  the  new  data. 

CALCULATIONS  AND  DATA  REDUCTION 

Figures  1  and  2  illustrate  the  method  used  to  compare  the  experimental 
results  with  theoretical  results  to  determine  the  empirical  continuum  in 
each  of  the  narrow  windows.  The  solid  curve  in  each  panel  is  based  on 
calculations  for  a  sample  with  the  indicated  equivalent  pressure  and 
absorber  thickness  to  match  an  experimental  sample.  Values  of  the  true 
transmittance  T'(calc)  for  the  solid  curve  were  calculated  according  to 
Equation  (8)  at  0.02  cm“l  intervals  by  summing  the  calculated 

contributions  by  all  of  the  lines  centered  between  2900  cm"l  and  4200 
cm"l.  Intensities,  widths  and  positions  of  the  lines  used  in  the 
computations  are  from  the  1982  version  of  the  AFGL'^  line  parameter 
tapes. 

The  upper  panel  of  Figure  1  covers  a  narrow  spectral  interval  on  the 
low-wavenumber  side  of  the  region  studied.  All  of  the  strong  absorption 

lines  are  centered  at  higher  wavenumbers,  and  most  of  the  absorption 

displayed  in  the  narrow  interval  between  3071.5  and  3072.5  cm’l  is  due 
to  the  extreme  wings  of  these  lines.  Only  a  very  small  fraction  of  the 
absorption  in  this  interval  is  due  to  the  two  weak  lines  centered  near 

3071.3  and  3072.8  cm“l.  The  short  curve  made  up  of  long  dashes  in  the 
upper  panel  of  Figure  1  is  from  the  experimental  spectrum  obtained  for  a 
laboratory  sample  with  the  same  parameters  as  those  used  for  the 
calculated  spectrum..  Because  of  the  finite  slitwidth  of  the  spectrometer 
used  to  obtain  the  data,  the  experimental  curve  is  smoothed  and  does  not 
show  much  of  the  structure  found  in  the  calculated  spectrum  of  T'(calc). 
In  order  to  make  a  more  direct  comparison  between  the  experimental  data 
and  theoretical  results,  we  computed  the  degraded  spectrum  of  T(calc), 
represented  by  the  curve  of  short  dashes,  by  convolving  the  calculated 
true  transmittance  spectrum  with  a  spectral  slitwidth  SW  =  0.6  cm"l  to 
match  the  spectral  slitwidth  of  the  spectrometer.  The  form  of  the 

triangular  slit  function  is  illustrated  in  the  lower  panel  of  Figure  1. 
The  appropriate  slitwidth,  which  varied  from  approximately  0.6  cm"^  near 
3000  cm“l  to  1  cm"l  near  4200  cm"^,  was  used  for  each  calculation  of 
T(calc). 
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Figure  1:  Spectral  plots  of  transmittance  in  two  narrow 

windows  near  3987.30  and  307Z.0  cm“^. 
Solid  curve,  calculated  true  transmittance; 
short'dashed  curve,  calculated  transmittance 
with  slit  function;  long-dashed  curve, 
experimental.  Upper  panel;  u  =  4.72  x 

10^2  molecules  cm"^,  p  =  0.0226  atm. 
Lower  panel:  u  =  1.37  x  10^2  molecules 

cm"2,  p  =  0.0226  atm. 
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Figure  2: 


Spectral  plots  of  transmittance  in  three 
narrow  windows  near  3682.78,  3717.32,  and 

3739.62  cm"l.  Solid  curve,  calculated  true 
transmittance;  short-dashed  curve,  calculated 
transmittance  with  slit  function;  long-dashed 
curve,  experimental.  u  =  3.66  x  10^^ 

molecules  cm"^,  p  =  0.0226  atm  for  all 
three  panels. 


The  curves  in  the  lower  panel  of  Figure  1  represent  a  narrow  window  on 
the  high  wavenumber  side  of  the  many  very  strong  lines  in  the  strong  band 
covered  in  this  study.  As  in  the  upper  panel,  much  of  the  calculated 
absorption  at  the  wavenumbers  of  maximum  transmittance  between  the  very 
weak  lines  is  due  to  the  extreme  wings  of  distant  lines.  Also,  as  in  all 
of  the  50  narrow  windows  investigated,  the  experimental  transmittance  near 
the  center  of  the  window  is  less  than  that  calculated  on  the  basis  of  the 
known  line  intensities  and  widths.  Within  each  window,  the  excess 
absorption  has  the  nature  of  a  continuum,  i.e.,  it  does  not  show  the 
structure  associated  with  lines.  Thus,  we  call  it  the  empirical  continuum 
because  it  provides  an  empirical  fit  between  the  theoretical  results  and 
the  experimental  data. 

The  value  of  the  normalized  empirical  coefficient,  for  self 

broadening  near  the  center  of  a  window,  is  determined  by  comparing  the 
experimental  value  T(exp)  near  the  wavenumber  of  maximum  transmittance 
with  the  corresponding  calculated  transmittance  T(calc),  which  includes 
convolution  of  the  slit  function. 

T(emp)  =  T(exp)  /  T(calc),  and  (14) 

®Cs  =  (-1/up)  (In  T  (emp))  (15) 

A  similar  method  is  used  to  determine  the  empirical  continuum 
coefficients  for  N2  broadening.  Samples  of  H2O  +  N2  used  for  this 

part  of  the  study  were  at  total  pressures  of  either  5  atm  or  10  atm  with 
the  H2O  partial  pressure  p  less  than  0.03  atm.  Thus,  the  self 
broadening  is  much  less  than  the  N2  broadening  for  these  samples,  and 
can  be  ignored  when  comparing  the  experimental  transmittance  to  the 
calculated  T(calc).  It  follows  that  the  corresponding  coefficient  for 
N2  broadening  is: 

=  (-1/upj^)  (In  T  (emp))  (16) 

This  method  of  comparing  experimental  results  with  the  calculated 
values  accounts  for  the  effects  of  the  finite  slitwidth  of  the 
spectrometer  and  for  the  very  weak  lines  that  occur  within  the  spectral 
interval  passed  by  the  slit.  Of  course,  errors  in  the  assumed  intensities 
of  these  very  weak  local  lines  lead  to  errors  in  the  coefficients 
determined  for  the  empirical  continuum.  However,  it  is  unlikely  that  such 
errors  are  serious  in  very  many  of  the  narrow  windows  investigated.  The 
intensitites  of  these  lines  have  been  determined  by  a  combination  of 
theory  and  experiment  prior  to  being  included  in  the  AFGL  listings,  and 
most  are  probably  accurate  to  within  a  few  percent.  In  most  of  the 
windows,  these  loca’  lines  contribute  only  a  small  fraction  of  the 
calculated  absorption;  therefore,  errors  in  their  assumed  intensities  are 
not  expected  to  cause  serious  errors  in  the  calculated  values  of 
transmi ttance. 

This  method  of  comparison  also  makes  it  possible  to  determine  the 
empirical  continuum  without  first  determining  the  true  values  of  the 
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absorption  coefficient  at  all  of  the  points  within  the  interval  passed  by 
the  spectral  slit.  The  validity  of  the  method  depends  on  the  assumption 
that  T(emp),  the  transmittance  of  the  empirical  continuum,  is  constant 
over  the  spectral  slitwidth.  The  results  indicate  that  this  assumption  is 
valid,  as  we  expect,  since  the  empirical  continuum  is  believed  due  to  the 
extreme  wings  of  the  lines. 

The  coefficients  ®CO  and  K®  are  normalized  to  one  atmosphere 
pressure  of  pure  H2O  for  self  broadening  and  the  same  pressure  of  N2 
for  N2  broadening.  Thus,  the  units  for  these  coefficients  include 
atm"l  along  with  the  molecules"^  cm"2.  Inclusion  of  the  atm"l 
indicates  that  -InT*  for  a  given  absorber  thickness  u  (molecules  cm"2) 
is  proportional  to  pressure.  In  accordance  with  the  discussion  in  Section 
1,  this  is  true  as  long  as  all  of  the  absorption  at  the  particular 
wavenumber  of  interest  is  due  to  the  wings  of  the  lines.  At  the 
wavenumbers  chosen  for  detailed  study,  this  is  essentially  true  for  p|\| 
less  than  one  atmosphere  and  for  the  H2O  partial  pressure  p  less  than 
0.2  atm.  Thus,  the  assumed  relationships  are  valid  for  applications  to 
the  earth's  atmosphere.  However,  at  higher  pressures  this  simple 
relationship  between  -InT'  and  pressure  may  not  hold  at  some  of  the 
wavenumbers  where  the  measurements  have  been  made.  This  is  true  because 
some  of  the  nearby  lines  are  broad  enough  that  their  half  widths  are  no 
longer  small  relative  to  the  distances  to  the  centers  of  the  lines.  (See 
Equation  (1).) 

Because  of  the  finite  spectral  width  of  the  sTit  used  in  the 
measurements,  the  spectrometer,  when  adjusted  to  the  center  of  one  of  the 
windows  of  interest,  may  receive  radiation  at  wavenumbers  that  include  the 
centers  of  some  of  the  nearby  lines.  Thus,  we  do  not  expect  the  measured 
-InT  to  be  exactly  proportional  to  up  (or  up^)  at  the  centers  of  these 
windows.  Comparing  the  experimental  results  as  we  did  with  the  calculated 
values  (including  the  convolved  slitwidth)  overcomes  the  problem  of 
measuring  the  coefficient  directly. 

The  method  used  also  enables  us  to  determine,  with  good  reliability, 
the  true  transmittance  T',  and  thus  the  corrected  value  of  K®,  that 
would  be  observed  with  zero  slitwidth.  This  corrected  coefficient  is 
given  for  self  broadening  by 

Kj  (corrected)  =  Ks(calc)  +  ®Cs  .  (17) 

A  corresponding  equation  applies  to  N2  broadening. 

As  a  check  against  possible  errors  in  our  method  of  reducing  the  data, 
we  investigated  the  absorption  in  some  of  the  windows  by  an  alternate 
method  that  depends  less  on  the  accuracy  of  the  assumed  intensities  of  the 
local  lines.  Values  of  the  empirical  continuum  coefficients  obtained  by 
this  method  agreed  well  with  the  others,  thus  adding  confidence  in  the 
accuracy  our  method  of  accounting  for  the  local  lines  and  the  finite 
slitwidth.  In  the  alternate  method,  we  obtained  experimental  spectra  for 
pure  H2O  samples  at  2  or  3  different  pressures,  and  plotted  the  quantity 
(-l/u)ln  T(exp)  vs  p.  In  accordance  with  Equation  (9)  and  the  discussion 
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following  the  equation,  the  slope  of  the  curve  through  the  plotted  points 
corresponds  to  the  experimental  coefficient  K(wing)  normalized  to  one  atm 
pressure.  From  this  experimental  value,  we  then  subtracted  the  calculated 
value  of  K(wing),  normalized  to  one  atm,  to  obtain  ®CO. 

Because  the  empirical  coefficients  are  derived  from  the  differences 
between  two  other  coefficients,  one  experimental  and  one  calculated,  the 
percentage  uncertainties  in  these  empirical  values  may  be  larger  than  for 
either  of  the  coefficients  from  which  they  are  derived.  Uncertainties  in 
the  experimental  values  vary  widely  from  one  portion  of  the  spectrum  to 
another.  The  largest  uncertainties  occur  in  the  wings  of  the  bands  where 
the  coefficients  are  so  small  that  even  the  largest  samples  produce  only  a 
small  amount  of  absorption.  The  uncertainty  is  also  larger  than  normal 
for  some  of  the  narrow  windows  in  the  stronger  part  of  the  band  where  the 
local  lines  make  a  larger  contribution,  or  where  the  observed 
transmittance  depends  strongly  on  the  slitwidth. 

The  estimated  uncertainty  in  the  empirical  coefficients  listed  in  the 
following  sections  increases  from  approximately  10%  or  15%  in  the  main 
part  of  the  band  to  as  much  as  50%  in  the  wings  of  the  band  where  the 
coefficients  are  small. 

SELF  BROADENING:  3000  -  4200  cm"! 

Figure  3  summarizes  the  results  for  self  broadening  at  the  centers  of 
approximately  50  narrow  windows  throughout  the  strong  vibration-rotation 
band.  The  calculated  coefficients  represented  by  the  circles  vary  in 
value  by  approximately  a  factor  of  1000  from  the  wings  of  the  band  to 
places  in  the  central  portion  of  the  band  within  a  few  cm"l  of  some  of 
the  very  strong  lines.  These  calculated  values  were  determined,  as 
explained  above,  by  summing  the  contributions  by  all  of  the  lines  while 
assuming  a  simple  Lorentz  line  shape. 

As  would  be  expected,  the  circles  do  not  fall  on  a  smooth  curve 
because  of  the  near-random  distribution  of  the  H2O  lines.  Some  of  the 
wavenumbers  may  be  within  2  or  3  cm"l  of  lines  that  contribute  much  of 
the  absorption,  whereas  others  may  not  be  within  several  cm"^  of  a  line 
that  contributes  significantly.  Furthermore,  the  intensities  of  lines 
within  an  interval  of  several  wavenumbers  may  vary  widely.  Note  the  two 
points  at  3775.2  and  3791.04  cm"^,  which  are  much  lower  than  the 
adjacent  points.  These  wavenumbers  occur  in  the  window  between  the  two 
branches  of  the  band. 

Corresponding  values  of  the  empirical  continuum  are  represented  in 
Figure  3  by  triangles  joined  by  lines.  In  the  wings  of  the  band,  the 
empirical  coefficients,  ®C^,  are  typically  from  3  to  6  times  as  large 
as  the  corresponding  calculated  coefficients.  This  ratio  is  smaller  for 
the  coefficients  near  the  very  strong  lines  around  3800  cm"^-  Each 
corrected  coefficient  is  determined  by  summing  the  corresponding 
calculated  coefficient  and  empirical  continuum  coefficient.  Thus,  the 
corrected  coefficient  represents  the  value  that  would  be  derived  from  a 
measuremnt  made  with  infinite  resolving  power  (zero  slitwidth).  In 
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Figure  3:  Spectral  plots  of  the  normalized  absorption  coefficients  from 

3000  to  4200  cm"l  for  self  broadening.  The  circles 
represent  monochromatic  transmittance  values  calculated  from 
line  parameters  and  the  Lorentzian  line  shape.  The  triangles 
represent  empirical  values  derived  from  the  ratio  of  the 
experimental  transmittance  values  to  the  monochromatic  values 
calculated  from  line  parameters  and  convolved  with  an 
instrument  slit  function  (see  Eqs.  14,  15).  The  corrected 
values  (solid  squares)  are  the  sums  of  the  empirical  and 
calculated  values  (see  Eq.  17).  Temperature,  296K. 


accordance  with  the  above  discussion,  such  a  measurement  must  be  made  for 
a  sample  at  a  pressure  sufficiently  low  that  the  half  widths  of  any  lines 
are  much  less  than  the  distance  from  the  point  of  measurement  to  the  line 
centers.  This  condition  is  necessary  for  -In  (T)/u  to  be  proportional  to 
pressure,  as  is  assumed  in  deriving  the  constants. 

The  left-hand  portion  of  Table  1  summarizes  the  three  sets  of 
coefficients  for  self  broadening  plotted  in  Figure  3:  the  calculated, 
corrected  and  empirical  continuum.  The  right-hand  portion  of  the  table 
includes  comparable  data  from  FASCODE,  a  widely  used  computer  program 
developed  at  Air  Force  Geophysics  Laboratory  for  line-by-line 
calculations.  These  data  are  discussed  at  the  end  of  this  sub-section. 

The  results  shown  here  should  not  be  mistakenly  interpreted  to  mean 
that  the  actual  absorption  coefficients  are  2  to  6  times  as  Targe  as 
calculated  values  everywhere  within  the  band.  This  is  definitely  not  the 
case.  Recall  that  the  wavenumbers  selected  for  the  measurements  were 
chosen  at  the  centers  of  narrow  windows  within  the  band.  Because  of  this 
method  of  selection,  the  absorption  coefficients  are  much  smaller  than 
they  are  only  a  few  cm"l  away  where  relatively  strong  lines  are 
centered.  Thus,  at  the  wavenumbers  studied,  a  large  fraction  of 
absorption  is  due  to  lines  centered  a  few  cm'l,  or  further,  away. 

Lines  centered  in  the  wings  of  the  band,  approximately  3000  -  3500 
cm"l  and  4000  -  4200  cm"^,  are  much  weaker  than  those  centered  in  the 
main  part  of  the  band.  A  calculated  coefficient  in  the  wing  of  the  band 
typically  contains  a  significant  portion  due  to  the  lines  centered  from 
1-10  cm"l  away,  but  they  also  contain  a  large  portion  due  to  more 
distant  lines.  Although  the  lines  centered  in  the  main  part  of  the  band 
are  very  far  from  the  narrow  windows  in  the  wings  of  the  band,  these  lines 
are  very  strong  and  therefore  make  a  sizeable  contribution  to  the 
absorption  in  the  narrow  windows  in  the  wings.  Thus,  a  significant  part 
of  a  coefficient  in  the  wings  may  be  due  to  lines  centered  more  than  100 
cm"l  away.  By  contrast,  the  narrow  windows  in  the  main  part  of  the  band 
are  much  nearer  the  very  strong  lines;  therefore,  only  a  small  fraction  of 
a  calculated  coefficient  in  this  region  is  due  to  lines  centered  more  than 
10  or  20  cm"l  away. 

We  attribute  the  empirical  continuum  to  a  "super  Lorentzian"  shape  of 
portions  of  the  wings  of  the  absorption  lines.  This  shape  can  be 
described  most  easily  in  terms  the  following  simple  equation: 

k  =  kL  X  (18) 

where  the  correction  factor  x  is  a  function  of  (v  -v  o)^  and  kL 
represents  the  coefficient  based  on  the  Lorentz  shape  given  by  Equation 
(1).  The  line  is  said  to  be  super-Lorentzian  when  x  >1,  and 
sub-Lorentzian  when  x  <1.  If  x  >1,  as  we  suggest,  for  the  wings  of 
the  lines,  it  must  be  less  than  unity  for  a  portion  of  the  line  near  the 
center,  in  accordance  with  Equation  (2).  Our  data  do  not  provide  any 
evidence  either  way  about  any  slight  deviations  from  the  Lorentzian  shape 
near  the  center  or  intermediate  wings.  However,  it  should  be  noted  that 
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only  a  small  fraction  of  the  integral  in  Equation  (2)  is  due  to  the  wings 
where  |  v  -  vq  1  is  greater  than  about  5  cm"l,  the  approximate  distance 
at  which  we  suspect  that  x  becomes  greater  than  unity.  Thus  x  would 
not  need  to  differ  much  from  unity  for  smaller  values  of  |v  -vqI  in 
order  for  the  integral  of  Equation  (2)  to  be  valid. 

It  could  be  informative  to  attempt  the  derivation  of  a  modified  shape 
for  the  lines  that  would  predict  the  corrected  values  of  the  continuum 
coefficient  listed  in  Table  1.  Such  a  task  is  beyond  the  scope  of  the 
present  study.  In  view  of  the  uncertainties  in  the  line  parameters  and 
the  experimental  data,  it  seems  unlikely  that  a  definitive  line  shape 
could  be  derived  to  fit  all  of  the  data.  Since  the  cause  of  the  deviation 
of  the  shape  is  not  understood,  it  might  not  be  realistic  to  assume  that 
all  of  the  lines  have  the  same  shape. 


Although  it  would  be  difficult,  if  not  impossible,  to  derive  a  single 
function  x  ('>"  ^o)  would  result  in  agreement  between  calculated  and 
observed  values,  certain  qualitative  characteristics  of  the  shape  seem 
consistent.  Results  of  the  present  study  and  previous  studies  of  this 
type  indicate  that  x  ^or  self-broadened  lines  becomes  greater  than  unity 
at  1  V  -  vq  I  equal  to  about  5  or  10  cm“^.  The  value  of  x  then 
probably  increases  with  increasing  j\>  -vqI  until  it  approaches  10  or 
more  at  (  v  -  vq  !  “  cm"l.  Recall  that  the  ratios  of  some  of  the 

corrected  coefficients  to  the  calculated  ones  in  the  wing  of  the  band 
below  3500  cm“l  are  between  10:1  and  5:1.  For  values  of  |v-vo  I  ^s 
large  as  100,  the  more  complex  line  shapes  discussed  in  references  5  and  6 
should  replace  the  overly  simple  Lorentz  shape  given  in  Equation  (1),  and 
empirical  changes  should  be  made  to  the  more  appropriate  shape. 

The  method  used  to  determine  the  empirical  continuum  coefficient  from 
data  obtained  with  spectral  slits  as  wide  as  1  cm'^  is  based  on  the 
assumption  that  the  continuum  is  constant  over  this  interval.  We  have 
checked  this  assumption  in  a  few  different  spectral  intervals,  and  the 
results  indicate  that  the  assumption  is  valid.  For  example,  we  consider 
the  interval  from  approximately  3680  to  3684  cm"^  shown  in  the  upper 
panel  of  Figure  2. 

At  3682,8  cm“^,  near  the  peak  transmittance,  the  experimental 

transmittance  is  0.400,  and  the  calculated  transmittance,  including  the 
slit  function,  is  0.710.  Dividing  the  first  of  these  values  oy  the 

latter,  we  find  T(emp),  the  transmittance  of  the  empirical  rontinuum,  to 
be  0.563,  Following  the  same  procedure,  we  find  T(emp)  to  be  0.570  at 
3679.5  and  0.556  at  3684.5  cm'l,  the  two  wavenumbers  where  the 
transmi ttances  are  at  minima.  In  addition,  we  found  the  average  T(emp) 
for  the  entire  interval  3679.5  -  3684.5  cm"^  to  be  0.546.  All  of  these 
values  for  T(emp)  agree  within  expected  experimental  error;  thus,  we 
conclude  that  this  quantity  is  essentially  constant  over  the  5  cm'l  wide 
i nterval . 


We  made  similar  checks  on  the  empirical  continuum  at  different  places 
within  several  other  intervals.  In  all  cases,  the  empirical  continuum 
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appeared  essentially  constant,  within  experimental  uncertainty,  over 
intervals  2-5  cm~l  wide. 


The  two  absorption  lines  centered  near  3679.5  and  3684.5  cm"l  are 
strong  enough  to  be  opaque  near  their  centers  for  the  sample  represented 
in  Figure  2.  However,  these  two  lines  contribute  only  a  small  fraction  of 
the  continuum  absorption  near  3682.8  cm"l.  Most  of  this  continuum 
arises  from  much  stronger  lines  centered  outside  of  the  interval  shown  in 
the  figure. 


Our  simple  suggested  model  of  a  super-Lorentzian  line  shape  would  lead 
to  an  empirical  continuum  that  changes  very  little  over  a  distance  of  a 


few  cm" 
exists 
greater 
that 
maximum 


1.  As 
because 

than  some  value,  approximately  5  to  10  cm"^.  It  is  likely 
X  increases  gradually  with  increasing  |  v-  \;ol  ;  it  may  reach  a 

and  start  decreasing,  but  the  changes  are  expected  to  be  gradual. 


explained  above,  we  believe  that  the  empirical  continuum 
X  in  Equation  (18)  is  greater  than  unity  for  |v  -  vqI 

10  cm-1. 


Thus,  the  sum  of  the  contributions  of  many  lines  would  change  gradually. 


Comparison  With  FASCODE  Scientists  at  AFGL  have  developed  a  very 
efficient  program  to  compute  absorption  coefficients  from  the  AFGL  line 
parameters^.  This  widely  used  program  called  FASCODE^  does  not  sum 

the  contributions  by  all  of  the  lines  in  the  inefficient  manner  employed 
for  the  calculations  discussed  above.  Instead,  it  calculates,  point  by 
point,  the  sum  of  the  contributions  of  the  nearby  lines,  those  within 
approximately  25  cm’^.  A  continuum  is  then  added  to  this  result  to 
account  for  the  more  distant  lines  and  for  empirically  determined 
deviations  from  the  simple  theoretical  line  shape. 


The  continuum  now  used  as  part  of  the  FASCODE  prograin  for  the  spectral 
region  covered  by  this  report  was  based  on  empirical  line  shapes  derived 
from  a  combination  of  theory  and  experimental  results  from  lower 
wavenumbers.  The  new  data  presented  in  this  report  provide  a  means  of 
checking  the  accuracy  of  the  FASCODE  continuum.  Note  that  the  FASCODE 
continuum  does  not  represent  the  same  physical  quantity  as  the  empirical 
continuum  represented  in  Figure  3  and  Table  1.  Our  empirical  continuum 
represents  the  difference  between  experimental  results  and  theoretical 
results  based  on  all  of  the  lines,  with  the  assumption  that  the  line 
shapes  are  given  by  Equation  (1).  On  the  other  hand,  the  FASCODE 

continuum  was  intentionally  made  to  include  the  contributions  due  to 
distant  (further  than  approximately  25cm~l)  lines  as  well  as  empirical 
corrections  for  the  super-Lorentzian  line  shapes  observed  at  lower 
wavenumbers. 


Figure  4  compares  the  corrected  coefficients,  based  on  our  self 
broadening  data,  with  the  corresponding  values  calculated  by  FASCODE. 
These  quantities  are  directly  comparable  because  they  relate  to  the  total 
coefficient,  which  includes  contributions  due  to  nearby  lines  and 
continuum.  The  corrected  coefficients  representing  our  data  in  Figure  4 
are  the  same  as  those  in  Figure  3.  We  note  that  in  most  cases  our 
corrected  values  are  greater  than  the  values  based  on  FASCODE  thus 
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indicating  that  the  FASCODE  continuum  should  be  increased  to  improve  the 
agreement  with  experimental  results. 


The  right-hand  portion  of  Table  1  gives  more  detailed  results  of  the 
FASCODE  calculations,  and  provides  a  point-by-point  comparison  with  our 
results.  Columns  3  and  5  contain  the  quantities  plotted  in  Figure  4. 
Column  6  shows  the  portions  of  each  coefficient  due  to  the  FASCODE 
continuum;  the  remainder  of  the  coefficient  in  column  5  is  due  to  the 
nearby  lines. 

In  almost  all  cases,  the  calculated  values  for  the  total  coefficient 
based  on  FASCODE  (column  5)  are  smaller  than  the  corresponding 
experimental  values  in  column  3  but  larger  than  our  calculated  values  in 
column  2.  The  latter  of  these  two  results  is  due  to  the  super-Lorentzian 
line  shape  assumed  by  the  authors  of  FASCODE  in  determining  their 
continuum.  Our  corrected  coefficients  indicate  that  the  lines  are  more 
super-Lorentzian  that  the  FASCODE  authors  assumed,  and  that  the  continuum 
should  be  increased.  The  values  suggested  for  the  continuum  are  listed  in 
the  last  column  of  Table  1.  Usage  of  these  suggested  values  instead  of 
those  now  used  would  bring  FASCODE  calculations  into  agreement  with  our 
experimental  results. 

N2  BROADENING:  3000  -  4100  cm"! 

Figure  5a  and  the  left-hand  portion  of  Table  2  summarize  the  results 
of  the  investigation  of  broadening  in  most  of  the  narrow  windows 
included  in  the  study  of  self  broadening.  As  in  Figure  3,  the  circles  in 
Figure  5a  represent  the  calculated  continuum  coefficients  based  on  the 
simple  Lorentz  shape.  The  normalized  halfwidths  of  ihe  lines  used  for 
N2  broadening  are  the  ones  in  the  AFGL  listing.  In  accordance  with  the 
self-broadening  coefficient  defined  by  Equation  (7),  these  half  widths  are 
one-fifth  as  large  as  the  corresponding  ones  for  self-broadening. 
Therefore,  the  calculated  continuum  coefficients  broadening  are 

smaller  by  the  same  ratio.  '  , 

The  overall  shapes  of  the  empirical  continuum  curves  for  self 
broadening  and  N2  broadening  are  similar.  However,  the  values  for  N2 
broadening  are  typically  less  than  one-tenth  the  corresponding  values  for 
self  broadening.  The  figures  and  Table  2  also  show  that  the  ratios  of 
the  empirical  continuum  coefficients  to  the  calculated  coefficients  are 
much  smaller  for  N2  broadening  than  for  self  broadening.  Also,  as  for 
the  self  broadening,  this  ratio  is  less  near  the  very  strong  lines  in  the 
main  part  of  the  band  than  it  is  in  the  wings. 

The  empirical  continuum  for  N2  broadening  can  also  be  explained  in 
terms  of  super-Lorentzian  lines.  The  factor  x  in  Equation  (18)  probably 
becomes  greater  than  unity  at  |  v  -  vg  'equal  to  approximately  10,  and 
continues  increasing  until  [v  -  vq  I  is  greater  than  100  cm“l. 
Measurements^"^  in  th*^  major  atmospheric  windows  indicate  that  x  for 
1N2  broadening  decreases  after  reaching  some  maximum  value  and  becomes 
less  than  unity  for  very  large  j  v-  jqI-  ^be  maximum  value  of  x  is 
apparently  much  greate*'  for  self-broadening  than  for  N2  broadening. 


Nj  BROADENING 
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Figure  5a:  Spectral  plots  of  the  normalized  absorption  coefficients  from 

3000  to  4200  cm"^  for  N2  broadening.  The  circles 
represent  monochromatic  transmittance  values  calculated  from 
line  parameters  and  the  Lorentzian  line  shape.  The  triangles 
represent  empirical  values  derived  from  the  ratio  of  the 
experimental  transmittance  values  to  the  monochromatic  values 
calculated  from  line  parameters  and  convolved  with  an 
instrument  slit  function  (see  Eqs.  14,  15).  The  corrected 
values  (solid  squares)  are  the  sums  of  the  empirical  and 
calculated  values  (see  Eq.  17).  Temperature,  296K.  (Compare 
with  Figure  3.) 
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Comparison  With  FASCODE.  Figure  5b  and  the  right-hand  portion  of 
Table  2  provide  comparisons  of  our  results  for  N2  broadening  with 
calculations  based  on  FASCODE.  Each  column  of  Table  2  corresponds  to  the 
same  column  for  self  broadening  in  Table  1. 

The  calculated  values  for  N2  broadening  based  on  FASCODE  generally 
agree  with  our  experimental  results  much  better  than  do  the  corresponding 
self  broadening  data  in  Table  1  and  Figure  5a.  At  many  of  the 
wavenumbers,  the  discrepancy  is  less  than  the  experimental  uncertainty. 


EMPIRICAL  CONTINUUM:  1300  -  2200  cm'l 

Like  the  3000  -  4200  cm’l  interval,  the  1300  -  2200  cm"l  region 
contains  many  very  strong  lines  that  make  up  a  vibration-rotation  band  of 
H2O  vapor.  This  region  is  bounded  on  both  sides  by  wide  windows,  and 
essentially  all  of  the  absorption  arises  from  lines  centered  within  the 
spectral  interval.  In  Reference  1,  we  reported  on  a  study  of  the 
continuum  absorption  in  several  narrow  windows  within  this  band.  The 
method  of  analysis  was  slightly  different  from  the  method  employed  in  the 
present  study  of  the  3000  -  4200  cm“l  region,  but  the  objectives  were 
the  same.  Because  of  the  strong  similarity  in  the  results  we  have 
included  three  figures  from  the  previous  report  for  comparison  with  the 
data  for  the  3000  -  4200  cm“l  region. 

Figure  6,  which  corresponds  closely  to  Figure  3  for  the 
higher-wavenumber  region,  shows  that  the  experimental  values  for  the 
self- broadening  continuum  coefficient  are  greater  than  the  calculated 
coefficients.  The  general  contour  of  the  empirical  continuum  curve 
follows  the  smooth  contour  of  the  absorption  band.  The  maxima  that  occur 
at  approximately  1530  cm"l  and  1659  cm"l  are  near  the  peaks  of  the  two 
branches  of  the  band,  and  the  minimum  near  1580  cm"l  is  close  to  the 
window  in  the  center  of  the  band. 

A  composite  of  several  curves  of  empirical  continuum  from  the  previous 
work  are  shown  in  Figure  7  for  both  self  broadening  and  N2  broadening. 
The  overall  shape  of  the  curves  for  N2  broadening  is  very  similar  to 
that  for  self  broadening,  but  the  largest  N2  broadening  coefficients  are 
less  than  one-tenth  as  large  as  the  corresponding  values  for  self 
broadening.  Both  of  these  results  are  very  consistent  with  the  findings 
in  the  3000  -  4200  cm"^  region. 

Figure  7  also  illustrates  an  obvious  strong  decrease  in  the  continuum 
coefficients  with  increasing  temperature,  another  result  that  is 
consistent  with  the  continuum  in  other  spectral  regions^"^’^. 


8.  D.  E.  Burch  and  R.  L.  Alt,  AFGL-TR-84-0128,  ADA  147391  Scientific 
Report  No.  1,  AFGL  Contract  No.  F19628-81-C-0118  (1984). 
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Figure  5b:  Corrected  absorption  coefficients  for  N2  broadening,  from 

Fig.  5a  and  Table  2,  compared  with  the  coefficients 
calculated  from  FASCODE  (see  Table  2),  including 
contributions  from  FASCODE  local  lines  and  continuum  (see 
reference  6).  Temperature,  296K.  (Compare  with  Figure  4.) 
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Spectral  plots  of  the  normalized  absorption  coefficients  from 
1300  to  1800  cm"i  for  pure  H2O.  The  empirical  continuum 
curve  is  drawn  through  the  +'s,  which  represent  the 
differences  between  the  experimental  values  (solid  squares) 
and  the  calculated  values  (circles)  bsed  on  the  line 
parameters.  Temperature,  308K.  (from  Ref.  1). 


empirical  continuum  at  various  temperatures, 
broadening;  and  N,  N2  broadening. 


A  quantitative  interpretation  of  the  obvious  decrease  in  the 
coefficients  with  increasing  temperature  is  complicated  by  the  shift  in 
the  intensities  of  the  lines.  As  the  temperature  increases,  the  lines  in 
the  wings  of  the  band  increase  in  intensity  at  the  expense  of  the  lines 
closer  to  the  band  center  (near  1595  cm"l).  This  shift  causes  a 
coefficient  in  the  wings  of  the  band  to  be  higher  at  high  temperatures 
than  it  would  be  if  the  lines  remained  constant  in  intensity  while 
changing  in  shape  and  width. 


SPECTRAL  BAND  CONTOURS  OF  LIQUID,  VAPOR,  AND  EMPIRICAL  CONTINUUM 

Some  research  workers  have  attributed  the  excess  absorption  by  H2O 
in  different  spectral  regions  to  clusters  of  HgO  molecules  that  take  on 
many  of  the  absorption  characteristics  of  liquid  water.  To  get  an 
indication  of  the  validity  of  this  hypothesis  we  have  compared  the  spectra 
of  the  empirical  continuum  with  those  of  liquid  water.  The  comparisons 
are  shown  in  Figures  8  for  the  3000  -  4200  cm"^  region,  and  in  Figure  9 
for  the  previous  work  in  the  1300  -  2200  cm"^  region.  Each  of  these 
figures  also  contains  a  curve  (V)  that  represents  the  smoothed  contours  of 
the  vibration-rotation  band  of  the  H2O  vapor. 

The  curve  labeled  (V)  in  Figure  8  was  obtained  by  summing  the 
intensities  of  all  the  absorption  lines  in  each  25  -  cm“l  wide 
interval.  The  value  to  be  plotted  for  a  given  interval  was  derived  by 
giving  full  weight  to  the  sum  for  that  interval  and  half  weight  to  each 
adjacent  interval.  These  weighted  values  were  then  plotted  to  form  the 
basis  for  the  curve.  A  similar  method  was  used  to  derive  the 
corresponding  curve  for  Figure  9.  These  two  curves  ’’abeled  (V)  show  the 
general  contour  of  the  band  without  the  excessive  structure  that  would 
exist  without  the  averaging. 

The  similarity  between  the  curves  for  vapor  and  empirical  continuum  in 
Figure  9  is  remarkable.  Both  curves  show  maxima  at  the  peaks  of  the  P- 
and  R-  branches  and  a  minimum  near  the  center  of  the  band.  In  the 
original  report^,  we  interpreted  this  similarity  as  evidence  that  the 
empirical  continuum  is  asociated  with  the  vibration-rotation  lines.  There 
is  little  similarity  between  the  empirical  continuum  curve  and  the  curve 
labeled  L,  which  represents  the  absorption  coefficient  of  liquid  water. 
Thus,  it  seems  unlikely  that  the  excess  absorption  represented  by  the 
empirical  continuum  is  due  to  liquid-like  clusters. 

The  similarity  between  the  vapor  curve  and  the  empirical  continuum 
curve  in  Figure  8  is  somewhat  less  than  it  is  in  Figure  9.  Both  curves  in 
Figure  8  are  highest  in  the  300  -  500  cm"l  wide  region  that  contains  the 
strongest  lines,  and  both  show  a  minimum  at  approximately  3701  cm'l  near 
the  band  center.  However,  the  empirical  continuum  also  contains  some 
structure  not  present  in  the  vapor  curve  between  3600  and  3800  cm"^,  a 
region  that  includes  several  very  strong  lines.  Structure  also  appears  in 
this  spectral  region  in  the  empirical  continuum  curve  for  N2  broadening 
in  Figure  5.  We  have  not  attempted  to  determine  quantitatively  the  cause 
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Figure  8:  Comparison  of  the  spectral  curves  between  3000  and  4200 

cm“l  of  the  empirical  continuum  for  self  broadening  (C), 
the  absorption  coefficient  of  liquid  water  (L),  and  the 
average  intensities  of  the  H2O  vapor  lines  (V).  The  arrows 
point  to  the  appropriate  scales  for  the  curves. 


of  this  structure,  which  is  uncharacteristic  of  empirical  continuum  curves 
in  other  spectral  regions 

The  similarities  in  the  shapes  of  curves  C  and  V  in  Figures  8  and  9 
suggest  that  the  excess  absorption  represented  by  the  empirical  continuum 
is  associated  directly  with  both  the  rotational  transitions  and  the 
vibrational  transition  involved  in  the  band.  It  seems  unlikely  that  the 
empirical  continuum  would  display  the  spectral  characteristics  it  does  if 
it  were  due  to  a  completely  different  process  that  involves  dimers  made  up 
of  two  bound  H2O  molecules  roatating  as  semi-rigid  rotators  or  vibrating 
in  a  mode  that  involves  relative  motion  between  the  two  molecules.  If  the 
excess  absorption  is  due  to  dimers,  as  been  suggested  by  other  research 
workers,  the  binding  is  apparently  such  that  each  H2O  molecule  forming  a 
dimer  is  still  free  to  vibrate  and  rotate  similar  to  the  way  it  does  in 
the  normal  monomer  state.  Significant  shifts  would  be  expected  in  the 
vibrational  and  rotational  frequencies  if  the  absorption  were  due  to 
dimers. 

The  very  significant  effect  of  the  N2  broadening  represented  in 
Figure  5a  and  7  suggests  further  that  the  excess  absorption  is  associated 
directly  with  the  vibration-rotation  band  rather  than  with  dimers  or 
clusters.  Continuum  absorption  by  dimers  or  clusters  is  expected  to  be 
essentially  independent  of  the  presence  of  nonabsorbing  N2.  Thus,  it 
seems  likely  that  much  of  the  excess  absorption  is  due  to  the  extreme 
wings  of  individual  absorption  lines  that  absorb  more  than  is  predicted  by 
any  of  the  widely  used  theoretical  line  shapes. 
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PREVIOUS  REPORTS  AND  PAPERS  UNDER  PRESENT  CONTRACT 


The  present  contract  has  dealt  almost  exclusively  with  the  infrared 
continuum  absorption  by  H2O  vapor.  Essentially  all  of  the  results  of 
the  experimental  investigation  are  included  in  this  report  and  a 
Scientific  Report  published  in  1984. 

The  results  of  some  of  the  work  was  also  presented  in  oral  papers  at 
the  1982  and  1983  Annual  Review  Conferences  at  AFGL. 

The  reference  and  the  abstract  of  the  Scientific  Report  are  as 
follows: 

D.  E.  Burch  and  R.  L.  Alt,  AFGL-TR-84-0128 
Scientific  Report  No.  1,  AFGL  Contract 
F19628-81-C-0118  (1984) 

New  laboratory  measurements  have  been  made  of  the  H2O  continuum 
absorption  at  several  wavenumbers  within  the  700-1200  cm“^  and 
2400-2800  cm"l  regions.  Samples  of  both  pure  H2O  and  H2O  +  N2 
have  been  studied  in  a  multiple-pass  absorption  cell  at  path  lengths 
up  to  more  than  800  meters.  Data  for  samples  of  pure  H2O  at  reduced 
temperature  (284  K)  indicate  that  the  self-broadening  coefficients 
increase  with  decreasing  temperatures  at  about  the  rate  predicted  by 
extrapolating  from  previous  data  on  samples  at  elevated  temperatures. 
Coefficients  for  N2  broadening  have  been  measured  in  both  windows, 
and  the  results  compared  with  previous  work.  The  results  have  been 
compared  with  the  data  base  for  the  well-known  transmission  computer 
code  LOWTRAN  6. 
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